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We discuss a method to select the velocities of ultra-cold atoms with a modified Fabry-Perot 
type of device made of two effective barriers and a well created, respectively, by blue and red 
detuned lasers. The laser parameters may be used to select the peak and width of the transmitted 
velocity window. In particular, lowering the central well provides a peak arbitrarily close to zero 
velocity having a minimum but finite width. The low-energy atomic scattering off this laser device 
is parameterized and approximate formulae are found to describe and explain its behaviour. 
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Velocity selection is a basic operation in quantum op- 
tics and atomic physics for a plethora of applications. 
There are mechanical (slotted disks) and non-mechanical 
(optical) techniques available, useful for different exper- 
imental circumstances, species, and energies. The large 
wavelengths achieved with laser cooling have made the 
traditional methods no longer effective because of the 
increasing importance of gravity and the quantum na- 
ture of translational motion. For example, the stan- 
dard classical-mechanical analysis of mechanical velocity- 
selection methods becomes invalid for small-time tem- 
poral slits, since they produce momentum spread in 
agreement with a time-energy uncertainty principle 0. 
Among the new methods, the velocity selection using 
Doppler sensitive stimulated Raman transitions J2| , and 
coherent population trapping into a dark state |3J, pro- 
vide selectivity in the "transverse direction" parallel to 
the lasers, and rely on specific internal level configura- 
tions. Fabry-Perot (FP) cavities have been also proposed 
to provide coherent velocity selection or trapping for lon- 
gitudinal motion, using detuned lasers perpendicular to 
the incident atoms Q or microwave cavities 0, |(| to im- 
plement the partially reflecting mirrors. The velocity se- 
lection in these cavities is produced by the filtering effect 
of resonance peaks in the transmission probability. The 
potential of FP cavities as trapping devices also stems 
from characteristic resonance features: high densities and 
large life times in the interaction region. 

The aim of this work is to discuss an improvement 
of these cavities, provide formulae to describe their be- 
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FIG. 1: m — mass( 23 Na); (a) Gaussian functions, d = 6pm, 
a — 2 fj,m; (b) square functions, d = 5 /im 
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FIG. 2: Transmittance versus velocity for different well 
depths. Vb = 300?i/s; V w = (dashed-dotted line), V w = 
150 h/s (thick dotted line), V w = 180.2 h/s (thick dashed line), 
V w = 180.25 h/s (solid line), V w = 180.5 h/s (thick solid line); 
other parameters in Fig. [Qi. The inset is a zoom of the 
lower-left corner. 

haviour, and study the fundamental limitations to lower 
the peak width and velocity of the transmitted wave 
packet. The basic idea is to add a well with control- 
lable depth between the two external barriers, see Fig. 
1. Effective barriers and well can be implemented with 
blue and red detuned lasers, respectively, which do not 
excite the impinging ground state atom and cause only a 
mechanical effect. The depth of the well and the barrier 
height can be varied with the intensities of the lasers. 
Making the well deeper, rather than wider, displaces 
the resonance peaks to lower energies without diminish- 
ing the inter-resonance spacing so it is the ideal way to 
achieve a sharp low-energy velocity selection. The reso- 
nance peak displacement with the well-depth can be seen 
in Fig. 2. The velocity shift is accompanied by a peak 
width reduction until a minimum, non-zero width is at- 
tained when the peak reaches zero velocity at a critical 
"threshold" depth. Beyond that depth the peak broad- 
ens, moves to higher velocities, and its maximum decays, 
as shown in the inset of Fig. 2. The effects of different 
depths are summarized in Fig. 3, which will be explained 
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FIG. 3: (a) Resonance velocity vr versus V w (symbols con- 
nected with dotted line). Filled symbols indicate the case 
|T(^' j r) 1 2 > 0.995, and empty symbols otherwise. The solid 
lines show the approximation of Eq. (J1J. Gaussian func- 
tions, see Fig. Vb = 3007i/s, a — 0.65 (diamonds), 
Vt — 500?i/s, a = 0.70 (squares); square functions, see Fig. 
GJ>: V b = 300ft/s, a = 0.79 (triangles up), V b = 500ft/s, 
a = 0.85 (triangles down), the circles indicate V w ,thres for the 
square model in the 2-pole approximation, (b) Velocity width 
Avr of the resonance versus V w ; meaning of symbols as in (a). 
The solid lines show the approximation of Eq. JKJ , the circles 
indicate (V w ,thres, Ai)R t thres) with the 2-pole approximation. 



next in more detail. 

We shall use both a realistic model based on three 
Gaussians, see Fig. 1, as well as a simplified version with 
two square barriers and a well. The scattering off the two 
potential models is very similar but the later enables us 
to obtain analytical exact results and approximate but 
physically illuminating expressions. Let us consider, for 
a single ultra-cold atom, the Hamiltonian 



H 



2m dx 2 



+ V b {x + d)-V w {x) + V b (x-d), (1) 



where V byW (x) = V bjW H(x), and II can take the forms 
■ 1 x r 1 if - d/2 < x < d/2 
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for the Gaussian and square models respectively. For 
simplicity we have set all Gaussians with the same width 
cr, and the square segments with the same length d. We 
assume that the atom impinges from the left and only ini- 
tial positive velocities are considered. We are interested 
in the transmission amplitude T and the "transmittance" 
\T\ 2 of the scattering solutions of H<f> v {x) = E v (f> v (x), 
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FIG. 4: Connection between peaks and poles in the square 
model, Vt = 400?i/s, other parameters in Fig. for the 
plotted parameter range, the 2-pole approximation and the 
exact result are indistinguishable; (a) resonance velocity vr 
(line with filled dots), resonance width Avr (line with un- 
filled diamonds) versus the distance ±A of the pole to the 
collision point (A < before the collision, A > after the col- 
lision); a 1-pole approximation with only ki is also plotted (vr 
(dashed line), Avr (dashed-dotted line)), (b) Motion of the 
two poles: V w = 199.884R/s (crosses), V w ,coii = 199.8987i/s 
(coinciding squares), V w ,thres ~ 199.9017i/s (small circles), 
V w = 199.927i/s (big circles). 



where E v = = Ti 2 k 2 /(2m). Both velocity, v, and 
wavenumber, k, will be used, the later being more ap- 
propriate for complex plane analysis and the former for 
presenting the physical results. For the square model, 



T(jfe) = ~Ae- Mk kk 2 k w (2) 
x | e ld ( fe - fe ») [ik b {k + k w )C+ (kk w - kl)S\ 2 

- e ld{k+k -\-ik b (k - k w )C+(kk w + A: fc 2 )^] 2 } 1 

where C = cosh(dk b ), S = sinh(dfc 6 ), k w = (k 2 + K 2 ) 1 / 2 , 
k b = [K 2 - k 2 ) 1 ' 2 , and K b = {2mV b ) 1 ' 2 /%, K w = 
(2mV w ) 1 / 2 /h. 

In both models, at some critical, "threshold" well- 
depths V Wt thres new bound states are formed, and for 
well-depths close to these thresholds, the description 
of the transmission peak is not as simple as for iso- 
lated Breit-Wigner (BW) resonances (see Fig. @J. At 
a depth V WtCO u slightly before threshold, a resonance- 
antiresonance pole collision occurs, and with further 
deepening, two "virtual" states in the complex momen- 
tum plane appear. Centering our attention on the first, 
lowest-energy resonance at zero well-depth, one can dis- 
tinguish when increasing the well-depth: first an ordinary 
resonance regime with a BW transmittance peak; second, 
an intermediate pre-bound state regime near threshold, 
in which the second pole cannot be ignored; and finally 
a bound-state regime. In Figure [2] transmittance curves 
corresponding to the different regimes are depicted. Even 
though the calculations can be made exactly, it is useful 
for applications and physically illuminating to describe 
these three stages in terms of approximate expressions 
and dominant dependences relating: poles of T in the 
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complex momentum plane; well-depth or other poten- 
tial parameters; and visible features such as position and 
width of the resonance. 

The first stage, dominated by a BW resonance pole k\ 
in the fourth quadrant of the momentum complex plane, 
is the most important one for velocity selection since it 
allows to diminish the resonance velocity and width by 
deepening the well, see Figs. |21and|2t- &i is accompa- 
nied by an antiresonance at ki — —k*. Since the bar- 
rier is symmetrical with respect to parity, the (Gamow) 
resonance states have well defined parity. Thus they ap- 
pear alternatively in one of the two eigenvalues of the 
2x25 matrix, for symmetrical, So, or antisymmetri- 
cal scattering, Si. The transmission amplitude is given 
by T = (S + Si)/2. We shall follow the motion of the 
first symmetrical resonance, the one that will become the 
ground state for deep enough wells, and assume that the 
first antisymmetrical pole of S\ is far from the origin so 
that a two pole approximation suffices. Retaining only 
two poles in the canonical pole expansion for cut-off po- 
tentials, 



So = - 



,-2ifcr (j; k*)(fc ^2) 

(k-h)(k-k 2 )' 



Si = e 



-2ikr 



(3) 



Here r is 3d/ 2 for the square model. For the Gaussian 
model, we could truncate the potential at a large r value 
and apply Eq. J3J). In any case the phase factor does not 
play any role to calculate the filtering function |T| 2 . 

Expressions for the two important poles can be ob- 
tained with the square model under some approxima- 
tions, as we shall see later on. In the "BW" regime the 
antiresonance ki may normally be ignored if the reso- 
nance is sharp (i.e. ki is close to the real axis) and far 
from the origin. A decrease in the well-depth displaces ki 
to the left and upwards, so that the transmittance curve 
decreases both its peak velocity and width. By inspec- 
tion of the S matrix, it is clear that in this regime the 
transmittance reaches the unitary limit \T\ 2 = 1 close to 
SRe(fci). We define Er = mv R /2 and vr as the energy 
and velocity of the transmittance maximum. In the BW 
regime Er ss h 2 $le(ki) 2 /2m. Simple parameterizations 
of this regime are provided by perturbation or semiclas- 
sical formulae. Let Erq and vrq be the real energy and 
velocity of the resonance peak "without well" (V w = 0). 
Then the energy of the resonance with non-zero well can 
be approximated within a perturbation theory for reso- 
nance functions Up to first order in V w , 



Er — E 



no 



aV w 



vr 



R0 



2aV w /m. (4) 



A semiclassical treatment |8j for opaque barriers gives 
a = 1, but keeping a as a fitting parameter the depen- 
dence of Eq. is valid even beyond very opaque barriers 
or very small depths, as can be seen in Fig. |3K- We define 
a velocity width Avr as the width of the transmittance 



peak at half height. In Fig. 03 the BW regime corre- 
sponds to the slow decrease with well-depth up to the 
abrupt, almost vertical increase associated with a bound 
state. A semiclassical estimate for the energy-width of 
the resonance is given by a well-frequency factor times 
the WKB probability to escape through a barrier from 
the well (see e.g. ||). Retaining dominant dependences 
in the opaque and shallow well limit, 



Av R = Avro exp(-/3T4,), 



(5) 



which, again, by keeping [3 as an effective fitting pa- 
rameter, describes the correct behaviour in the whole 
BW regime, until well-depths very near the intermedi- 
ate threshold region, see Fig. Eh- 

Near the threshold depth the velocity and width of the 
peak are affected more and more by the nearby antireso- 
nance, hi — —k\. This intermediate regime is extremely 
narrow, with respect to variations of V w , compared to the 
BW and bound state ones (see Fig. 0}. Nevertheless, its 
analysis is worthwhile since it establishes the ultimate 
physical lower limit of the peak velocity and width us- 
ing a FP filtering device. At a critical "collision" depth 
V w ,coll both poles meet at -iK co u, k coU > 0, on the neg- 
ative imaginary axis (see e.g. the "square" in Fig. 
in one dimensional scattering, as for s-wave scattering, 
the collision is not at the origin because bound states are 
not degenerate). Note that, in spite of their zero real 
part, the velocity peak is not at v — 0. As the well be- 
comes more profound the two poles move in opposite di- 
rections, now along the imaginary axis as "virtual" poles 
until the upper one arrives at the origin at the threshold 
depth V w jhres, with the lower pole at —inthres (see e.g. 
"circles" in Fig. @J|. The motion of the two poles just 
before the collision and even beyond threshold is well de- 
scribed by expanding the denominator of So in powers 
of (V w — Vw^coii) and retaining the first term. This gives 
fci,2 = -iKcoll ± h(Vw ~ V WtCO u) 1/2 i with 7 real. A con- 
sequence is that, at threshold, ki = 0, ki = —inthres ~ 
—2in co ii. The threshold is a singular, abnormal point in 
which the transmission peak reaches the origin, T(0) = 1 
(T(0) = for any other well-depth). Moreover, from the 



2-polc approximation of S, \T\ 2 = K 2 hres /(k 



^thres 



Thus, the width at half height, considering only positive 
momenta, reaches its minimum value. In wavenumber 
units it is just Kthres, arid the maximum, \T\ 2 = 1, occurs 
at k = Kthres- Approximate expressions for Kthres may 
be obtained from Eq. J5J. We have to find zeros of the 
denominator of T. Let x = We assume Vb,V w > 
and xCl. Neglecting 0(x 3 ) we arrive at the quadratic 
equation a2X 2 + 2iaix ~ a o = with ao — cot (^K w ^j — 
K w coth (dK b ) /K b , ai = K w / (2K b sinh 2 (dK b )), and 
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K 2 (dK w + sin (dK w )) 
AKl sin 2 (dK w /2) 



4 



K w (coth (dK b ) (cosh 2 (dK b ) - 3) + dK b ) 
+ 2K b sinh 2 (dK b ) ' 

The two solutions are given by X1/2 = — icti/ofc ± 
y/a a 2 - a\/a 2 - At V w , thres = ^ K. 2 wthre j2m, a zero 
of the denominator of T is at k = 0, so a = 0. Note 
that if K b 3> 1, K W: thres ~ n/d. The other pole is at 
—inthres — — 2iK b a\ / a 2 and determines the minimal ve- 
locity width Av R j hres = TiKthres/m (see Fig. EJ). 

Finally, with further well deepening, the upper pole 
crosses the real axis and becomes a bound state. As 
predicted by Eq. |(3J, the transmittance peak broadens 
dramatically and moves to higher positive velocities; also 
the peak maximum becomes smaller than one, so this 
regime is no longer useful for velocity filtering, see Figs. 
121 and El 

As an application example we shall compute the trans- 
mitted velocity distributions resulting from the FP fil- 
tering of an an atomic wavepacket prepared as a Bosc- 
Einstein condensate 0, 0] in a trap. The trap is moved 
with a certain velocity with respect to the laboratory 
frame and turned off suddenly at t = 0. The con- 
densate expands until the nonlinear interaction between 
the atoms can be neglected and encounters the FP cav- 
ity. (Alternatively the triple potential can be moved 
with the trap at rest.) First we calculate numerically 
the ground state 4>o (normalized to 1) in the reference 
frame where the trap is at rest. A harmonic trap is as- 
sumed with frequency oj x in x direction and lu vz in y 
and z directions. Using a one-dimensional approxima- 
tion of the Gross-Pitaevskii equation the Hamiltonian is 

H $) = l£& + ^ x2 + 2hNauj yz \4>(x)\ 2 \ where N is 
the number of atoms in the condensate and a the scat- 
tering length. We take a = 2.93 x 10~ 9 m. Then we 
change to the lab frame where the trap moves with ve- 
locity vq. The ground state ipo in this reference frame 
is tf^oix) — e lx ~^~'4>Q{x). At t = the trap is turned 
off being at position xtrap and the velocity selection 
potentials are switched on, i.e., the time-evolution is 
given by the Gaussian version of the Hamiltonian of 
Eq. plus a term V representing the decaying non- 
linear effect due to free expansion in y and z directions, 
V = 2hNauj yz |V>(x)| 2 /(l + w 2 ^ 2 ). Fig. shows the 
momentum distribution at t = (ground state). At 
t = 0.8 s the non-linearity has practically vanished, so 
the momentum distribution stays stable until the veloc- 
ity selection. The filtered distributions at t — 8 s for 
several resonance velocities obtained with different well- 
depths are also shown. 

We have in summary proposed an improvement of 



Fabry-Perot cavities to select the velocity of ultra-cold 
atoms using a well between the partially reflecting mir- 
rors, and have provided simple formulae to explain and 
describe their behaviour and the minimal velocity peak 
(zero) and width (non-zero) that can be achieved. 
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FIG. 5: |V(M)| 2 for v = 0.0336 cm/s, xtrap = -600 ^m; 
V b = 300ft/s, lo x = 5/s, Lo yz = 100/s,iV = 5 x 10 4 ; t = 
(dotted line); t = 0.8s: solid line; t = 8s: V w = 140ft/s (thick 
dotted line), V w = 150?i/s (dashed line), V w = 160?i/s (thick 
solid line); the circles mark the resonance velocities vr. 
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